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Abstract 

As the luminosity of an accreting black hole drops to a few percent of Eddington, the spectrum switches from the 
familiar soft state to a hard state that is well-described by a distended and tenuous advection-dominated accretion 
flow (ADAF). An ADAF is a poor radiator, and the ion temperature can approach 10 12 K near the center, although 
the electrons are cooler, with their temperature typically capped at ~ 10 9 ~ n K. The foundational papers predicted 
that the large thermal energy in an ADAF would drive strong winds and jets, as later observed and also confirmed 
in computer simulations. Of chief interest, however, is the accreting gas that races inward. It carries the bulk of the 
accretion energy as stored thermal energy, which vanishes without a trace as the gas passes through the hole's event 
horizon. One thus expects black holes in the ADAF regime to be unusually faint. Indeed, this is confirmed by a 
comparison of accreting stellar-mass black holes and neutron stars, which reside in very similar transient X-ray binary 
systems. The black holes are on average observed to be fainter by a factor of ~ 100 — 1000. The natural explanation 
is that a neutron star must radiate the advected thermal energy from its surface, whereas a black hole can hide the 
energy behind its event horizon. The case for an event horizon in Sagittarius A* , which is immune to caveats on jet 
outflows and is furthermore independent of the ADAF model, is especially compelling. These two lines of evidence 
for event horizons are impervious to counterarguments that invoke strong gravity or exotic stars. 
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1. Introduction 

The foundations of our present understanding of 
advection-dominated accretion were laid out in a 
series of papers by Narayan & Yi (1994, 1995a, b, 
hereafter NY94, NY95a, NY95b), Abramowicz et al. 
(1995) and Chen et al. (1995), although some ideas 
were anticipated much earlier by Ichimaru (1977). 
The specific abbreviation ADAF, which stands for 
advection-dominated accretion flow, was introduced 
by Lasot£{I](1996) and has become standard in the 
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field. We review here the application of ADAFs to 
accreting black holes (BHs), and in §4 we also con- 
sider their application to accreting neutron stars. 

2. Accretion Regimes 

The energy equation per unit volume of an accre- 
tion flow can be written compactly in the form 
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where p is the density, T is the temperature, S is 
the entropy per unit mass, t is time, v is the flow 
velocity, and q+ and g_ are the heating and cooling 
rates per unit volume. This equation states that the 
rate at which the entropy per unit volume of the 
gas increases is equal to the heating rate minus the 
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cooling rate. Since any entropy stored in the gas is 
advected with the flow, the left-hand side of equation 
(1) may be viewed as the effective advective cooling 
rate (fodv Equation (1) can then be rearranged to 
give 

9+ = <7- +<7adv, (2) 

which states that the heat energy released by viscous 
dissipation is partially lost by radiative cooling and 
partially by advective cooling. 

The standard thin accretion disk model (Shakura 
& Sunyaev 1973; Novikov & Thome 1973; Frank, 
King & Raine 2002) corresponds to the case when 
the accreting gas is radiatively efficient, so that we 
have 

Thin Disk: q_ > q adv , L ~ O.lMc 2 . (3) 

Since the gas cools efficiently, the sound speed is 
much less than the local Keplerian speed and 
the disk is geometrically thin. Also, the disk radiates 
about a tenth of the rest mass energy of the accreting 
gas, the precise fraction depending on the BH spin 
(see Shapiro & Teukolsky 1983). 

An ADAF corresponds to the opposite regime. 
Here the gas is radiatively inefficient and the accre- 
tion flow is underluminous. Thus, an ADAF is de- 
fined by the condition 

ADAF: g ad v>g-, i«0.1Mc 2 . (4) 

Some papers in the literature define an ADAF as a 
flow that corresponds exactly to a self-similar solu- 
tion described in NY94 (see §3.1). This is a needless 
restriction. In our opinion, it is more fruitful to cm- 
ploy the general definition of an ADAF as given in 
equation (g]0]- 

There are two distinct regimes of advection- 
dominated accretion. The first, the one that we 
focus on in this article is when the accreting gas is 
very tenuous and has a long cooling time (NY94; 
NY95b; Abramowicz et al. 1995). This regime is 
sometimes referred to as a RIAF — a "radiatively 
inefficient accretion flow" — and is defined by the 
condition 

ADAF/RIAF : t coo i > t acc . (5) 



2 Actually, at mass accretion rates for which the flow is close 
to the boundary between the thin disk and ADAF solutions, 
we have (by continuity) g at j v Q- , and the flow is only 
marginally radiatively inefficient (see Fig. 4). Nevertheless, 
even here, the two solutions are very distinct from each other. 



Here t coo \ is the cooling time of the gas and t acc is 
the accretion time. 

The second regime of advection-dominated accre- 
tion is when the particles in the gas have no trou- 
ble cooling, but the scattering optical depth of the 
accretion flow is so large that the radiation is un- 
able to diffuse out of the system. This radiation- 
trapped regime was briefly discussed by Begelman 
(1979) and was developed in detail by Abramowicz 
et al. (1988) in their "slim disk" model. The defining 
condition for this regime of ADAFs is 

ADAF/Slim Disk : t dis > t acc , (6) 

where tdie is the diffusion time for photons. 

The present review is devoted exclusively to 
the ADAF/RIAF form of accretion (see Narayan, 
Mahadevan & Quataert 1998b; Kato, Fukue & 
Mineshige 1998; Lasota 1999a,b; Quataert 2001; 
Narayan 2002, 2005; Igumenshchev 2004; Done, 
Gierlinsky & Kubota 2007 for reviews emphasizing 
various aspects of ADAFs) . We will henceforth drop 
the modifier RIAF and refer to these flows simply 
as ADAFs. 

As explained above, an ADAF is by definition very 
different from the standard thin accretion disk. Cor- 
respondingly, it is characterized by very distinct ob- 
servational signatures. Both ADAFs and thin disks 
have well-known counterparts in nature. 

3. Properties of ADAFs 

3.1. Geometry and Kinematics 

In an ADAF, since most of the energy released 
by viscous dissipation is retained in the gas, the 
pressure is large, and so is the sound speed: c s ~ 
vk, where the Keplerian velocity vk(R) is equal 
to c/\/2r with r = R/Rs being the radius in 
Schwarzschild units, Rs — 2GM/c 2 . The large 
pressure has several immediate consequences. First, 
the accretion flow becomes geometrically thick, 
with a vertical height H of order the radius R 
(an ADAF may be viewed as the viscous rotating 
analog of spherical Bondi accretion). Second, the 
flow has considerable pressure-support in the radial 
direction, so the angular velocity becomes sub- 
Keplerian. Third, the radial velocity of the gas is 
relatively large: v ~ avx (H/R) 2 ~ ac/r 1 / 2 , where 
a ~ 0.1 — 0.3 (see below) is the standard dimen- 
sionless viscosity parameter (Shakura & Sunyaev 
1973). Fourth, the large radial velocity leads to a 
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short accretion time: t acc = R/v ~ ts/a, where 
t s = (2GM/R 3 ) 1 / 2 is the free-fall time. Finally, the 
large velocity and large scale height cause the gas 
density to be very low, and so the cooling time is 
very long and the medium is optically thin. 

The above properties are nicely illustrated in the 
self-similar ADAF solution derived by NY94, in 
which all quantities in the accretion flow behave as 
power-laws in radius. (A similar solution was ob- 
tained by Spruit et al. 1987 in a different context.) 
NY94 obtained a general solution for arbitrary vis- 
cosity parameter a, adiabatic index 7, and advec- 
tion parameteiE] / = <Zadv/9+- In the limit a 2 « 1 
(a good approximation) and / — > 1 (radiatively 
very inefficient flow), the solution simplifies to 
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The numerical coefficients on the right correspond 
to 7 = 1.5, a reasonable choice (see Quataert & 
Narayan 1999). 

The great virtue of the above self-similar solution 
is that it is analytic and provides an easy and trans- 
parent way of understanding all the key properties 
of an ADAF. Its biggest deficiency is that it is scale- 
free, which means that it is inappropriate near the 
inner or outer boundary of the flow. Therefore, for 
detailed work, one must use global solutions of the 
ADAF equations that satisfy appropriate boundary 
conditions (e.g., Abramowicz et al. 1996; Narayan, 
Kato & Honma 1997c; Chen, Abramowicz & La- 
sota 1997; Manmoto, Mineshige & Kusunose 1997; 
Popham & Gammie 1998; Manmoto 2000). Calcu- 
lating global solutions is somewhat involved and one 
may wish to employ some short-cuts (e.g., Yuan, 
Ma & Narayan 2008). On the other hand, even the 
exact global solutions are somewhat limited since 
they solve a set of height-integrated equations. For 
greater realism, one might wish to work directly with 
numerical simulations (e.g., Goldston, Quataert & 
Igumenshchev 2005; Noble et al. 2007). 



3.2. Thermal Properties 

The ADAF solution is gas pressure dominated. 
Since c s ~ Vk, this means the gas temperature is 
nearly virial. Under normal conditions, gas at such 
a high temperature will radiate copiously, especially 
at small radii where the temperature can approach 
10 12 K. Thus, in order to have an ADAF, the accret- 
ing gas generally has to be a two-temperature plasma 
(at least at small radii), with electron temperature 
T e much less than the ion temperature Tj (NY95b) . 
(The only way to avoid this condition is by having 
an extremely low accretion rate below about 10~ 6 
of the Eddington rate.) Typical ADAF models have 
the two temperatures scaling as 



Ti ~ 10 12 K/r, T e - Min (T h 10 9_li K) 



(10) 



3 By this definition, / = corresponds to a fully cooling- 
dominated (no advection) flow and / = 1 corresponds to a 
fully advection-dominated (no radiative cooling) flow. 



In order for gas in an ADAF to be two- 
temperature, there must be weak coupling between 
electrons and ions. Models generally assume that 
the coupling occurs via Coulomb collisions, which 
is inefficient at the densities under consideration. 
Begelman & Chiueh (1988) investigated whether 
plasma instabilities might enhance the coupling 
and drive the plasma rapidly to a single temper- 
ature; this would be problematic for the ADAF 
solution. However, they were unable to identify a 
clear mechanism. For normal mass accretion rates, 
the electrons in a two-temperature ADAF will have 
a thermal energy distribution (but not necessarily 
the ions, see Mahadevan & Quataert 1997). 

Early work on two-temperature ADAFs assumed 
that viscous heating acts primarily on the ions; for 
instance, NY95b took the ratio of electron heating 
to total heating, S, to be zero, while Esin, McClin- 
tock & Narayan (1997) assumed S ~ m e /m p ~ 
10~ 3 . However, an ADAF does not require 5 to be 
this small. Because of a degeneracy in model pa- 
rameters (Quataert & Narayan 1999), it is possible 
to have a viable ADAF model with larger values of 
8, provided the mass loss parameter s, defined in 
§3.6, is adjusted. More recent ADAF models (e.g., 
Yuan, Quataert & Narayan 2003) typically assume 
8 ~ 0.3-0.5. 

Various attempts have been made to estimate the 
value of 5 from first principles by considering the ef- 
fects of reconnection (Bisnovatyi-Kogan & Lovelace 
1997; Quataert & Gruzinov 1999) or MHD turbu- 
lence (Quataert 1998; Blackman 1999; Medvedev 
2000) . These studies do not agree on a single value 
of 8, but generally suggest that 8 is likely to be much 



3 



larger than 1CT 3 . Recently, Sharma et al. (2007) 
considered heating by the dissipation of pressure 
anisotropy and showed that S ~ (T e /T i ) 1 / 2 /3. 

It should be noted that, even if electrons and ions 
receive equal amounts of the dissipated energy (i.e., 
S = 0.5), the plasma can still be two-temperature. 
This is because a large part of the heating in an 
ADAF is by compression (since the density increases 
inward). Once kT > m e c 2 , which is the case for r < 
a few hundred, the electrons become relativistic and 
have an adiabatic index j e ~ 4/3, whereas the ions 
continue to be non- relativistic with 7, ~ 5/3. Since 
adiabatic heating by compression causes the tem- 
perature to scale as T oc p 7_1 , the electrons heat 
up only as T e <~ p 1 / 3 whereas the ions heat up as 
Tj - p 2 / 3 . Thus, even in the limiting case of S ~ 0.5, 
ADAFs naturally become two-temperature at small 
radii. For instance, in the ADAF model of Sagit- 
tarius A* (Sgr A*) proposed by Yuan et al. (2003), 
the authors obtain T e <~ O.ITi close to the BH even 
though they assumed 6 = 0.55. 

An important property of the ADAF solution is 
that it is thermally stable (NY95b; Wu & Li 1996; 
Kato et al. 1997; Wu 1997). The demonstration of 
this property was a crucial advance. Two decades 
earlier, in a seminal paper, Shapiro, Lightman & 
Eardley (1976), and after them Rees et al. (1982; 
'ion tori'), introduced the idea of a two-temperature 
accretion flow and derived a hot two-temperature 
accretion flow solution, the SLE solution. However, 
that solution turned out to be thermally unstable 
(Piran 1978). Until the development of the ADAF 
solution, and the recognition that it is different from 
the SLE solution, no stable, hot, optically thin solu- 
tion was available to model the many accretion sys- 
tems whose spectra (especially in the hard state, see 
below) demand such a flow. 

Chen et al. (1995) and Yuan (2003) have explored 
the relationships among the ADAF/RIAF, SLE, 
ADAF/Slim Disk and Thin Disk solutions. 

3.3. Models and Spectra 

The ADAF is a full accretion solution which incor- 
porates consistent dynamics, thermal balance, radi- 
ation physics, etc. Therefore, the radial profiles of 
all gas properties can be calculated self-consistently 
once we know the values of certain parameters: BH 
mass M, accretion rate M, viscosity parameter a, 
pressure parameter j3 = -P gas /(-Pgas + -Pmag), adia- 
batic index 7 (usually <~ 1.5), viscous heating pa- 



rameter 6, advection parameter /. 

Actually, apart from the system-specific param- 
eters M and M (which may be estimated through 
observation), most of the other parameters are con- 
strained. Under the nearly collisionless conditions 
expected in an ADAF, the viscosity parameter is 
moderately enhanced relative to a collisional gas 
(Sharma et al. 2006). In the case of dwarf no- 
vae in the hot state, Smak (1999) estimates a ~ 
0.2, while numerical simulations of the magneto- 
rotational instability give a < 0.1 (Hawley, Gammic 
& Balbus 1996). Thus we expect a ~ 0.1 - 0.3 
for an ADAF (typical values used in models are 
0.2 — 0.3). Numerical simulations further suggest 
that magnetic fields are generally subthcrmal, with 
-Pmag/-fgas ~ 0.1 (e.g., Hawley et al. 1996), so we 
expect (3 <~ 0.9. By calculating the energy loss via 
radiation from the hot accretion flow (synchrotron, 
bremsstrahlung, Compton scattering), / can be 
obtained self-consistently (e.g., NY95b; Narayan, 
Barret & McClintock 1997a; Esin et al. 1997; Yuan 
et al. 2003). Thus, we have only one poorly con- 
strained parameter: 5. The recent work of Sharma 
et al. (2007) provides a serviceable prescription even 
for this parameter. At large radii, where the plasma 
is effectively one-temperature, their formula gives 
S <~ 0.3, while in the energetically important inner 
region, where the plasma is two-temperature, they 
find S ~ 0.01 — 0.1, depending on model details. 
The final two parameters are discussed later: the 
transition radius R tT , §3.5, and the wind parameter 
s, §3.6. 

Since the earliest days of X-ray astronomy, it has 
been clear that BH binaries (BHBs) have a number 
of distinct spectral states (see Zdziarski & Gierlin- 
ski 2004; McClintock & Remillard 2006; Done et al. 
2007) . The most notable among these are the lumi- 
nous high soft state, or thermal state, the slightly 
less luminous low hard state, and the very under- 
luminous quiescent state. 

The thermal state is well described by the thin 
disk model, and a multi-color disk (MCD) black- 
body model (e.g., diskbb, Mitsuda et al. 1984; 
ezdiskbb, Zimmerman et al. 2005; both available in 
XSPEC, Arnaud et al. 1996) has been successfully 
used for years to model the X-ray spectra of sources 
in this state. Recently, fully relativistic versions of 
the MCD model for arbitrary BH spin (kerrbb, Li et 
al. 2005; bhspec, Davis & Hubeny 2006) have been 
developed, based on the relativistic thin disk model 
of Novikov & Thorne (1973). These models provide 
excellent fits to the X-ray spectra of BHBs in the 
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Fig. 1. Combined TTM (2-20 keV), HEXE (20-200 keV), 
and OSSE (50-600 keV) spectrum of GRO J0422+32 ob- 
served between 1992 August 29 and September 2. The solid 
line is an ADAF model which reproduces the overall lumi- 
nosity of the source, the power-law slope below the peak and 
the shape of the high-energy cutoff. The model has a mass 
accretion rate ~ O.lMEdd- (From Esin et al. 1998) 



thermal state (e.g., McClintock et al. 2006; Davis, 
Done & Blaes 2006). 

Whereas a satisfactory theoretical model, viz., the 
thin disk model, was established early on for the 
thermal state, the hard state was for many years a 
mystery. Figure 1 shows the spectrum of a typical 
BHB, GRO J0422+32, in the hard state (Esin et al. 
1998). The observations indicate that the accreting 
gas is very hot, kT > 100 keV. The gas must also be 
optically thin, since optically thick blackbody emis- 
sion at a temperature of 100 keV would correspond 
to a luminosity L = <jT a A > 10 46 ergs -1 for any 
reasonable estimate of the radiating area A. 

The most natural explanation of the emission in 
J0422 and other BHBs in the hard state is that it 
is produced by thermal Comptonization. Until the 
ADAF model was established, no accretion model 
could reproduce such a spectrum. (The SLE solution 
could, but it was unstable.) Indeed, astronomers 
were reduced to using empirical Comptonization 
models (Sunyaev & Titarchuk 1980) in which they 
postulated a Comptonizing cloud with some arbi- 
trary geometry and parameterized the cloud with an 
adjustable temperature and an optical depth (e.g., 
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Fig. 2. The spectrum of the BHB source V404 Cyg in qui- 
escence. The three dots on the left represent optical fluxes 
(Narayan et al. 1996), the error box in the X-ray band cor- 
responds to 2cr limits from ASCA (Narayan et al. 1997a), 
and the upper limit in the EUV is derived from the absence 
of a Hell A4686 line. The solid line corresponds to an ADAF 
model with M ~ 0.005MEdd- The dashed line is a thin disk 
model whose M has been adjusted to fit the optical data. 
This model fits poorly in the X-ray band and is inconsistent 
with the EUV limit. (From Narayan et al. 1997a) 

Zdziarski et al. 1996, 1998; Gierlinski et al. 1997). 

The situation changed with the recognition of the 
ADAF solution. This model turned out to have the 
precise properties — density, electron temperature, 
stability — needed to explain the hard state. The 
solid line in Fig. 1 shows an ADAF model of J0422 
(Esin et al. 1998) in which the accretion rate has 
been adjusted to fit the spectrum; the required rate 
is about a tenth of the Eddington mass accretion 
rate MEdd, where MEdd = £Edd/(0.1c 2 ), i.e., it is 
the mass accretion rate at which a disk with radia- 
tive efficiency 0.1 would radiate at the Eddington 
luminosity. It is gratifying that both the temper- 
ature (which determines the position of the peak) 
and the Compton y-parameter (which determines 
the power-law slope below the peak) are reproduced 
well. Esin et al. (1997, 1998, 2001) present models 
of other BHBs in the hard state. 

A typical accreting BH observed in the hard state 
has a luminosity on the order of a few percent of 
Eddington. At much lower luminosities, we have the 
quiescent state, where the spectrum becomes notice- 
ably different. Figures 2 and 3 show observations 
of a quiescent BHB, V404 Cyg, and a quiescent su- 
permassive BH, the Galactic Center source Sgr A*. 
Although these spectra look very different from the 
one shown in Fig. 1, the ADAF model is able to fit 
these and other observations of quiescent systems 
(Narayan, McClintock & Yi 1996; Narayan et al. 
1997a; Yuan et al. 2003). All it requires is a lower 
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Fig. 3. The quiescent spectrum of Sgr A*. The radio data are 
from Falcke et al. (1998; open circles) and Zhao et al. (2003; 
filled circles), the IR data are from Serabyn et al. (1997) 
and Hornstein et al. (2002), and the two "bow-ties" in the 
X-ray band correspond to the quiescent (lower) and flaring 
(higher) data from Baganoff et al. (2001, 2003). The solid 
line is an ADAF model of Sgr A* in the quiescent state. The 
mass accretion rate is M ~ IO^Me^ near the BH. (From 
Yuan et al. 2003) 

value of M, as appropriate for the lower luminos- 
ity. The qualitative features of the spectrum, e.g., a 
softening of the X-ray power-law index (see Corbel, 
Tomsick & Kaaret 2006), follow naturally. However, 
a caveat is in order: In some cases, the X-ray emis- 
sion in quiescence may be from a jet lauched from 
the ADAF (§3.6), rather than from the ADAF itself 
(e.g., Yuan & Cui 2005). 



3.4. Radiative Efficiency 

As expressed in equation (j4|), the defining prop- 
erty of an ADAF is that it is radiatively inefficient, 



ADAF 



r\ = L/Mc 2 < 0.1. 



(11) 



Calculations show that the ADAF solution is pos- 
sible only for low mass accretion rates. Specifically 
only when L < O.lLEdd is the gas density low enough 
to permit a two-temperature plasma (NY95b; 
Narayan 1996; Esin et al. 1997). Near the critical 
luminosity i cr ;t or critical mass accretion rate M cr i t 
at which the ADAF solution first becomes viable, 
the radiative efficiency is, by continuity not very 
different from that of a thin disk: rj ~ 0.1. However, 
with decreasing M, the efficiency decreases. Very 
roughly, we estimate (Fig. 4) 
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Fig. 4. Radiative efficiency r\ of an accretion flow around a 
BH, as a function of the Eddington-scaled mass accretion 
rate (Hopkins, Narayan & Hernquist 2006b). The horizontal 
segment between 0.01 and 0.1 of the Eddington accretion rate 
shows a transition regime in which a part of the accretion 
flow is an ADAF, but the radiative efficiency is still large. 
It might correspond to the Intermediate State, and perhaps 
the upper end of the Hard State (§3.5, Fig. 5). Although this 
plot is based on calculations shown in Fig. 11 of NY95b and 
Fig. 13 of Esin et al. (1997), it is still very approximate. 

r) ~ 0.1 (M/0.01M E dd) , M < 0.01M Edd . (12) 

(We note that for the different prescription for i5 
used by Sharma et al. 2006, n falls rapidly only for 
M < 10- 4 M Edd .) 

The extreme inefficiency of an ADAF at very low 
accretion rates is critical for understanding the pe- 
culiar properties of accreting BHs in the quiescent 
state (as first discussed by Narayan et al. 1996 for 
BHBs and Narayan et al. 1995 for supermassive 
BHs). We do not discuss this topic here, but point 
the reader to other reviews for details (e.g., Quataert 
2001; Narayan 2002, 2005). The quiescent state has 
also played a major role in our efforts to test for the 
presence of an event horizon in BHBs. This is a key 
topic of §4. 

3.5. Spectral Regimes 

Based on the properties of the ADAF solution 
discussed above, Narayan (1996) proposed a simple 
model for understanding the spectral states of ac- 
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Fig. 5. Configuration of the accretion flow in different spec- 
tral states, shown schematically as a function of the mass 
accretion rate M (based on Esin et al. 1997). The ADAF 
is represented by the hatched ellipses, with the intensity of 
hatching indicating the density of the hot gas. The horizontal 
lines represent a standard thin disk. The lowest panel shows 
the quiescent state (QS), which corresponds to a very low 
mass accretion rate (say < 10 — Mfidd); a weak jet (§3.6), 
a low radiative efficiency (Fig. 4), and a large transition ra- 
dius (Fig. 6). The second panel from the bottom shows the 
hard state (HS), where the mass accretion rate is higher 
(~ 10~ 3 - 10 _1 - 5 M Bdd ) but still below the critical rate 
M cr i t , the jet is stronger, the radiative efficiency is some- 
what larger, and the transition radius is smaller. The second 
panel from the top shows the intermediate state (IS), where 
M ~ 10 -1 - 5 - 10 _1 M Edd < M crit , the jet is even stronger, 
the radiative efficiency is high ~ 0.1, and the transition ra- 
dius is fairly close to the ISCO. Finally, the top panel shows 
the thermal state (TS), where there is no ADAF, the thin 
disk extends down to the ISCO, r\ ~ 0.1, and there is no jet. 
Esin et al. (1997) included a tentative proposal for the so- 
called very high state (see also Done et al. 2007), now called 
the steep power-law state (SPL, McClintock & Remillard 
2006), but we do not include this still- mysterious state here. 

crcting BHs. This picture was developed in detail 
by Esin et al. (1997). The basic idea is illustrated in 
Fig. 5. 

The key parameter that determines the spectral 
state of an accreting BH is M. When M > M a it, 
only the thin disk solution is available, and so the 
accretion flow is in the form of a thin disk all the way 
down to the innermost stable circular orbit (ISCO). 
The system is then in the thermal state, and its 
spectrum is well-described by the MCD model. 

Once M falls below M cr ;t, both the thin disk and 
ADAF solutions become viable (at least at small 
radii). Now accretion continues as a thin disk at radii 



larger than a transition radius, R > R tI , but the 
flow switches to an ADAF at smaller radii. When 
the transition from a pure thin disk (thermal state) 
to a disk-plus- ADAF configuration first occurs, i.e., 
when M is just below M C rit, the ADAF is very small 
in size and we have a more-or-less radiatively effi- 
cient flow, with r\ <~ 0.1. This corresponds to the 
so-called intermediate state. Then, at a somewhat 
lower M, the ADAF expands a bit and rj is mod- 
estly lower and we have the classic hard state. Fi- 
nally, when M is much lower than M cr it , the ADAF 
becomes much larger, R tr ~ 10 3 — 10 i?s> with r/ <C 
0.1 (Narayan et al. 1996, 1997a; Menou, Narayan & 
Lasota 1999b). This is the quiescent state. 

The above paradigm has proved durable (see 
Zdziarski & Gierlinski 2004; Done et al. 2007). In 
particular, considerable evidence has accumulated 
that the thin disk retreats from the innermost stable 
circular orbit (ISCO) to a large radius in the qui- 
escent state. The evidence is strongest in transient 
BHBs and CVs. The spectra of quiescent BHBs 
show absolutely no sign of any soft blackbody- 
like X-ray emission from a thin disk at small radii 
(Narayan et al. 1996, 1997a; McClintock, Narayan 
& Rybicki 2004). Timing properties also indicate 
a large a large size for the ADAF (Hyncs et al. 
2003; Shahbaz et al. 2005). In additional, theoret- 
ical arguments indicate that the thermal-viscous 
disk instability, which causes the transient behavior 
in these systems, is incompatible with observations 
unless the disk is severely truncated in the quies- 
cent state (Lasota, Narayan & Yi 1996b; Hameury, 
Lasota & Dubus 1999; Lasota 2001, 2008; Dubus, 
Hameury & Lasota 2001; Yungelson et al. 2006). 
There are also indications from the time delay be- 
tween the optical and X-ray outbursts in the BHB 
GRO J1655-40 (Orosz et al. 1997; Hameury et al. 
1997) that the cool disk is truncated at a large 
radius in the quiescent state. 

In the case of supermassivc BHs in the quiescent 
state, there is no feature in the spectrum that might 
be associated with a thin disk, suggesting that there 
is no disk at all; examples are Sgr A* (Narayan et 
al. 1998a; Narayan 2002) and M87 (Di Matteo et 
al. 2000, 2003). Intermediate luminosity AGN do 
exhibit optical emission from a disk, but the "big 
blue bump" is much less pronounced than in high- 
luminosity AGN (Ho 1999); this suggests that the 
disk is truncated at a radius <~ 10 — 100Rs and the 
interior is filled with an ADAF (Gammie, Narayan & 
Blandford 1999; Quataert et al. 1999). Incidentally, 
the thermal- viscous disk instability does not appear 
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to operate in AGN disks (Menou & Quataert 2001; 
Hameury, Lasota & Viallet 2007). 

In the more luminous hard state, again, there is 
considerable spectral evidence that the disk is trun- 
cated at a transition radius outside the ISCO, and 
that the inside is filled with an ADAF-like hot flow. 
The most spectacular example is the BHB XTE 
J1118+480, for which observations carried out in 
the hard state had unprecedented spectral coverage. 
The observations are fit well with an ADAF model, 
with a transition radius at ~ 50i?s (Esinetal. 2001). 
The model even fits the complicated timing behavior 
of the source (Yuan, Cui & Narayan 2005). It is hard 
to imagine that the same data could be explained 
with any model in which a cool disk extends down to 
the ISCO. Done et al. (2007) review spectral obser- 
vations of a number of other BHBs where again the 
data require a truncated disk. Nemmen et al. (2006) 
show that both the spectrum and the double-peaked 
Balmer line profile of the LINER source NGC 1097 
are consistent with a disk truncated at a few hun- 
dred Rs- 

In an interesting study of Cyg X-l, Gilfanov, Chu- 
razov & Revnivtsev (1999; see Cui et al. 1999; Done 
et al. 2007; for discussions of other sources) found 
that, as the characteristic frequency in the variabil- 
ity spectrum of the source increases, the power-law 
tail in the spectral energy distribution steepens and 
the amplitude of the reflection component in the 
spectrum increases. This is exactly what one expects 
when the transition radius between the outer cool 
disk and the inner hot ADAF varies. With decreas- 
ing transition radius, (i) the noise frequency (which 
is likely related in some fashion to the Keplerian fre- 
quency at the transition radius) should increase, (ii) 
the hot medium should be cooled more effectively by 
soft photons from the disk, giving a steeper power- 
law tail, and (iii) the solid angle subtended by the 
cool disk at the ADAF should increase, and there 
should be a larger reflection component. Zdziarski, 
Lubinski & Smith (1999) have noted that a corre- 
lation between spectral slope and reflection is com- 
monly seen in both BHBs and AGN. 

A few BHBs in the hard state have been found to 
show a soft blackbody-like component in their spec- 
tra (Balucinska-Church et al. 1995; Di Salvo et al. 
2001; Miller et al. 2006a,b; Ramadevi k Seetha 2007; 
Rykoff et al. 2007) . This could be interpreted as ev- 
idence that the thin disk extends all the way down 
to the isccEl Some authors have noted that it is 



4 However, Di Salvo et al. (2001) estimate from their ob- 



possible for a thin disk to evaporate to an ADAF at 
a relatively large radius and for the hot gas to then 
re-condense into a thin disk at small radii (Rozan- 
ska & Czerny 2000; Liu et al. 2007; Mayer & Pringle 
20070 Such models might explain the occurrence 
of a soft spectral component in some hard state 
sources. However, we note that the soft component 
typically has only 10% of the total observed lumi- 
nosity. It is hard to understand how a radiatively 
efficient thin disk located at the ISCO could be such 
a minor contributor to the emitted radiation. 

This difficulty is highlighted in the work of 
D'Angelo et al. (2008). They point out that, in any 
model of the hard state, there will be considerable 
interaction between the hot gas which produces the 
hard X-rays and cool gas that may be present in a 
conventional disk. The interaction will occur via X- 
ray irradiation (unless one has large outward beam- 
ing, which seems unlikely, e.g., Narayan & McClin- 
tock 2005) as well as particle bombardment. Using 
a prototype model for ion bombardment (Deufel 
et al. 2002; Spruit & Deufel 2002; Dullemond & 
Spruit 2005), D'Angelo et al. (2008) show that a 
weak soft component in the X-ray spectrum arises 
quite naturally when the thin disk is truncated at 
~ 15 — 20 <J 6 1. The model fits the observations 
surprisingly well; in fact, the same model would 
predict a much stronger soft component, and would 
strongly disagree with the observations, if the cool 
disk were to extend down to the ISCO. 

One of the most vexing problems in the theory 
of ADAFs is that we do not have a robust method 
of estimating the location of the transition radius 
i?tr- From the earliest studies (e.g., NY95b; Meyer 
& Meyer-Hofmeister 1994), it has been plausibly ar- 
gued that i? tr increases with decreasing M. How- 
ever, reliable predictions of the exact dependence 
have proved difficult, although a number of studies 
have come up with semi-quantitative results (Meyer 
& Meyer-Hofmeister 1994; Honma 1996; Rozanska & 
Czerny 2000; Meyer, Liu & Meyer-Hofmeister 2000; 
Spruit & Deufel 2002; Mayer & Pringle 2007). Yuan 



servations of Cyg X-l that the disk inner edge is located at 
tens of Rg, not at the ISCO 

5 Models of this kind may give hard spectra without a sig- 
nificant decrease in radiative efficiency. This would be con- 
sistent with the transition region between mass accretion 
rates of 0.01 and 0.1 of Eddington shown in Fig. 4, where 
we have suggested an ADAF with high radiative efficiency 
may be present. 

6 This is a reasonable location for the transition radius in 
the hard state (Esin et al. 1998; Done et al. 2007) 
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Fig. 6. Estimated location of the dimensionless transition 
radius Rt T /Rs, plotted as a function of the Eddington-s- 
caled accretion luminosity L/L^ddi f° r a sample of BHBs 
and low-luminosity AGN (Yuan & Narayan 2004). The indi- 
vidual estimates of R tI are obtained by fitting spectral ob- 
servations and are very uncertain. Nevertheless, there seems 
to be a trend of increasing transition radius with decreasing 
luminosity. 



& Narayan (2004) tried to use observations to de- 
duce the run of Rt T with accretion luminosity. Fig- 
ure 6 shows their results. 

Another issue that has been discussed recently 
is hysteresis in the transition between the thermal 
state and the hard state. Specifically, with increas- 
ing M, the hard state survives up to fairly high lu- 
minosities ~ 0.1 — 0.2LEdd, whereas with decreas- 
ing M, the transition from the thermal state to the 
hard state occurs at a much lower luminosity ~ 
0.02 - 0.05L E dd (Miyamoto et al. 1995; Maccarone 
& Coppi 2003; Zdziarski et al. 2004; Done et al. 
2007). Meyer- Hofmeister, Liu & Meyer (2005) have 
provided a plausible explanation. Their disk evapo- 
ration model, coupled with Compton-cooling of the 
hot electrons in the ADAF, naturally produces a 
hysteresis in the location of the transition radius. 

The brightest systems in the hard state (L > O.lLEdd) 
are difficult to model with the standard ADAF 
model. A variant of the ADAF solution — a natural 
extension of the model — called the luminous hot 
flow (LHAF; Yuan 2001, 2003; Yuan & Zdziarski 
2004; Yuan et al. 2007; see also Machida, Nakamura 
& Matsumoto 2006), looks promising for modeling 
these sources. 



3.6. Outflows and Jets 

NY94, NY95a discovered an unexpected property 
of the ADAF solution: the accreting gas in an ADAF 
has a positive Bernoulli varametei 7 l which is de- 
fined as the sum of the kinetic energy, potential en- 
ergy and enthalpy, 

Be = v 2 /2-GM/R + w. (13) 

A positive Bernoulli constant means that the gas 
is not bound to the BH. (This is not surprising, 
since the gas is not losing energy through radiation.) 
Therefore, the above authors suggested that ADAFs 
should be associated with strong winds and ietcf^l 
(see also Meier 2001). 

Strong outflows have been seen in numeri- 
cal simulations of ADAFs. The first indications 
came from 2D and 3D hydrodynamic simulations 
(Stone, Pringle & Begelman 1999; Igumenshchev 
& Abramowicz 2000; Igumenshchev, Abramowicz 
& Narayan 2000), but it was soon confirmed in 
MHD simulations as well (Stone & Pringle 2001; 
Hawley & Balbus 2002; Igumenshchev, Narayan 
& Abramowicz 2003; Pen, Matzner & Wong 2003; 
Machida, Nakamura & Matsumoto 2004; McKinney 
& Gammie 2004; Igumenshchev 2004). 

Apart from producing a gas-dominated, large- 
scale outflow, MHD simulations of ADAFs also have 
a second distinct outflow component along the axis 
in the form of a collimated, Poynting-dominated, 
relativistic jet (McKinney 2005, 2006). The orig- 
inal suggestion of NY94 that ADAFs would have 
outflows and jets has thus been confirmed by these 
simulations. Nevertheless, some authors have dis- 
puted a connection between a positive Bernoulli 
parameter and an outflow since it is possible to 
come up with explicit models that have a positive 



7 Technically, it should be called the Bernoulli "function" 
since the quantity varies with radius, but we use Bernoulli 
"parameter" since that was the term used in NY94 

8 There has been a tendency in the literature to ignore these 
early papers and to credit Blandford & Begelman (1999) for 
establishing the connection between ADAFs and outflows. 
We note that the abstract of the first paper on ADAFs, 
NY94, states: "Further, the Bernoulli parameter is positive, 
implying that advection-dominatcd flows are susceptible to 
producing outflows... We suggest that advection-dominatcd 
accretion may provide an explanation for... the widespread 
occurrence of outflows and jets in accreting systems." In 
the second paper, NY95a, the title itself mentions outflows, 
"Advection-Dominated Accretion: Self-Similarity and Bipo- 
lar Outflows," and an entire paragraph of the abstract is 
devoted to the topic. 
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Bernoulli parameter but no outflow (Paczynski 
1998; Abramowicz et al. 2000). 

Observational evidence for the association of non- 
thermal relativistic jets with ADAFs has accumu- 
lated in recent years with the discovery of radio emis- 
sion in virtually every BHB in the hard state (Cor- 
bel et al. 2000; Fender 2001; Fender, Belloni & Gallo 
2004; Fender & Belloni 2004). The radio emission 
is generally too bright to be produced by thermal 
electrons in the accretion flow. It is therefore very 
likely to come from nonthermal electrons in a jet; 
in fact, radio VLBI imaging has revealed a resolved 
jet in a few sources. Thus, it is now observationally 
well-established that the hard state/ ADAF is asso- 
ciated with relativistic jets. According to the discus- 
sion in §3.5, the quiescent state also has an ADAF 
and should have a (weaker) jet. Indeed, radio emis- 
sion has been seen from two quiescent systems, V404 
Cyg (Hjellming et al. 2000; Gallo, Fender & Hynes 
2005) and A0620-00 (Gallo et al. 2006), confirming 
this expectation. 

How much of the X-ray emission in an ADAF sys- 
tem comes from the accretion flow and how much 
from the jet? A strong correlation has been seen 
between radio and X-ray luminosity in the hard 
state and quiescent state (Corbel et al. 2003; Gallo, 
Fender & Pooley 2003). At first sight this might sug- 
gest that the X-ray must also come from the jet (e.g., 
Falckc, Kording & Markoff 2004). However, since 
the jet flows out of the ADAF and is thus highly 
coupled to it, any model in which the X-ray emis- 
sion is from the ADAF and radio is from the jet is 
equally compatible with the observations. In recent 
times, several authors have come down in favor of 
an ADAF origin for the X-rays (e.g., Heinz & Sun- 
yaev 2003; Merloni, Heinz & Di Matteo 2003; Heinz 
2004; Heinz et al. 2005; Yuan et al. 2005). 

Additional evidence on this issue comes from hard 
state spectra indicating thermal emission (e.g., see 
Fig. 1), which is very different from the nonthermal 
power-law emission one expects from a jet. Markoff, 
Falckc & Fender (2001), among others, have argued 
that the X-ray emission is due to synchrotron emis- 
sion from a carefully tuned distrubution of nonther- 
mal electrons. However, Zdziarski et al. (2003) and 
Zdziarski & Gierlinski (2004) showed that the high- 
energy cutoff in the spectra of BHBs in the hard 
state, which arises naturally in a thermal ADAF- 
like model, is very difficult to reproduce in a non- 
thermal synchrotron model. A variant of the jet 
model assigns the X-ray emission to thermal elec- 
trons in the "base of the jet" (Markoff, Nowak & 



Wilm 2005), but the discussion then becomes se- 
mantic. The base of the jet is surely embedded in the 
underlying ADAF and it is not clear that one gains 
anything by simply relabeling the radiation from 
the ADAF as jet emission. The issue is discussed in 
greater depth in Narayan (2005), who advocates a 
'jet- ADAF' model (as described in Yuan et al. 2005; 
Malzac, Merloni & Fabian 2004), in which the high 
energy X-ray emission in the hard state comes from 
the ADAF and the low-energy radio (and infrared) 
emission comes from a jet. 

The situation is less clear-cut in the quiescent 
state, where nonthermal emission from the jet might 
dominate even in X-rays (Yuan & Cui 2005; Wu, 
Yuan & Cao 2007). It is interesting, however, that 
the most quiescent system we know, Sgr A* , shows 
no sign of a jet. Radio images with a resolution of 
~ 15i?s appear to be jet-free (Shen et al. 2005), 
while the quiescent X-ray emission is spatially re- 
solved and appears to be from thermal gas near the 
Bondi radius at <~ 10 5 Rs, not from a jet (Baganoff 
et al. 2001, 2003). 

Turning our attention now to the extended wind 
that is predicted from an ADAF, and seen in nu- 
merical simulations, one consequence of the wind is 
that, at each radius, a fraction of the accreting gas 
is lost from the system. Thus, the accretion rate it- 
self varies with radius. This is usually parameter- 
ized with an index s such that M(r) oc r s . Opera- 
tionally, this means that, in the self-similar regime, 
the gas density varies with radius as p oc r ~ 3 ^ 2+s 
rather than as r~ 3//2 in the NY94 model. Unfortu- 
nately, apart from the rather general constraint that 
s should lie between (no mass loss) and 1 (the 
limit of a CDAF, §3.7), there is no good theoreti- 
cal estimate of the value of s. We have to resort to 
numerical simulations or observations. The former 
tends to give somewhat larger values, e.g., s ~ 0.7 
(Pen et al. 2003), while the one case where the lat- 
ter approach has been tried, Sgr A* , yields a smaller 
value, s ~ 0.3 (Yuan et al. 2003). 

3.7. Convection 

Another general property of ADAFs was high- 
lighted by NY94 (see also Begelman & Meier 1982): 
ADAFs have entropy increasing inward and should 
be convectively unstable by the Schwarzschild cri- 
terion. Numerical hydrodynamic simulations of 
ADAFs confirmed the presence of strong convec- 
tion (Igumenshchev et al. 2000; Narayan, Igumen- 
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shchev & Abramowicz 2000) and led to the devel- 
opment of an analytical self-similar model called 
the convection-dominated accretion flow (CDAF; 
Narayan et al. 2000; Quataert & Gruzinov 2000). 

The CDAF model employs a simplified one- 
dimensional treatment of the accretion flow, in 
which all fluxes are assumed to be in the radial di- 
rection. The resulting density varies as p oc r -1 / 2 . 
In the language of §3.6, a CDAF corresponds to 
s = 1. However, technically, an idealized CDAF has 
no mass outflow, only an outward flow of energy by 
convection; the energy is assumed to flow out into a 
surrounding medium. In practice, the flow is never 
perfectly one-dimensional; 2D effects intrude and 
one expects substantial winds and mass loss as well. 
A real ADAF thus involves a complicated interplay 
between convection and winds. 

Magnetic fields introduce further complexity, and 
there has been some discussion in the literature on 
whether MHD ADAFs do or do not have real con- 
vection (Machida, Matsumoto & Mineshigc 2001; 
Balbus & Hawley 2002; Narayan et al. 2002b; Pen et 
al. 2003; Igumenshchev 2004). The question has no 
practical significance, however, since there is general 
agreement that ADAFs have both unstable entropy 
gradients and strong magnetic stresses. 

In the particular case of a spherical, non-rotating, 
MHD flow (the magnetic analog of the Bondi prob- 
lem), some complications are absent and the prob- 
lem becomes relatively clean. Here one finds that 
convection and magnetic fields strongly influence 
the accretion physics, and the resulting flow is very 
different from the standard Bondi solution (Igumen- 
shchev & Narayan 2002; Igumenshchev 2004, 2006). 
This result is likely to have significant implications 
for astrophysics; for instance, isolated neutron stars 
accreting from the interstellar medium will be very 
much dimmer than one might expect based on the 
Bondi solution (Perna et al. 2003). 

3.8. ADAFs Around Supermassive Black Holes 

The ADAF solution is essentially independent of 
the mass M of the central BH. That is, if length 
and time are scaled by M and the accretion rate 
is scaled by the Eddington rate (also proportional 
to M), then the same solution is valid for any BH 
mass. Therefore, any successful application of the 
ADAF model to a BHB system has immediate conse- 
quences for supermassive BHs (SMBHs) in an equiv- 
alent state, and vice versa. This close connection 



between ADAFs in BHBs and ADAFs in AGN was 
highlighted in Narayan (1996); it is also empirically 
obvious from Fig. 6, which combines observations 
and models of BHBs and AGN. 

The first SMBH to be modeled as an ADAF 
was the Galactic Center source Sgr A* (Narayan 
et al. 1995). This was soon followed by Fabian & 
Rees (1995), who suggested that quiescent nuclei in 
nearby giant ellipticals (e.g., M87) must be accret- 
ing via ADAFs, and by Lasota et al. (1996a), who 
argued that Low Ionization Nuclear Emission-line 
Region sources (LINERs, e.g., NGC 4258) and low- 
luminosity active galactic nuclei (LLAGN) must 
have ADAFs. Both suggestions have turned out to 
be correct (e.g., Reynolds et al. 1996; Mahadevan 
1997; Gammie et al. 1999; Quataert et al. 1999; Di 
Matteo et al. 2000, 2003; Loewenstein et al. 2001; 
Ulvcstad & Ho 2001; Nemmen et al. 2006). 

Other authors have suggested that all of the fol- 
lowing systems have ADAFs: FRI sources (Baum, 
Zirbel & O'Dea 1995; Reynolds et al. 1996; Begel- 
man & Celotti 2004), BL Lac sources (Maraschi & 
Tavecchio 2003), X-ray Bright Optically Normal 
Galaxies (XBONGs; Yuan & Narayan 2004), and 
even some Seyferts (Chiang & Blaes 2003). All of 
these sources are relatively low-luminosity AGN, 
where an ADAF is likely to be present (Figs. 5, 6). 

As discussed in §3.6, a feature of the ADAF model 
is the presence of outflows and jets. This implies that 
ADAF systems should generally be radio-loud. This 
is certainly the case for the FRI and BL Lac sources 
mentioned above. Ho (2002; see also Nagar et al. 
2000; Falcke et al. 2000) presents in Fig. 5b of his pa- 
per a very interesting plot of radio loudness, defined 
as the ratio of the flux at 6 cm to the flux in the opti- 
cal B band versus Eddington-scaled luminosity, for 
a sample of galactic nuclei. Virtually every source 
in the plot with an Eddington-ratio below about 
0.01 is radio loud, which is perfectly consistent with 
our expectation that all such sources should have 
ADAFs. Moreover, the degree of radio loudness in- 
creases with decreasing Eddington ratio, again con- 
sistent with our expectation that the ADAF should 
become more and more dominant with decreasing 
M. The majority of sources in Ho's plot that have 
L/L,Edd > 0.1 are radio-quiet, as we would expect 
from Fig. 5 if these systems have cool disks. How- 
ever, a small minority of these bright sources do show 
powerful jets (they are generally FRII sources). The 
exact nature of the accretion in these sources is un- 
clear. A more complete plot, with updated data, can 
be found in Sikora, Stawarz & Lasota (2007). 
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The jet and the extended wind from an ADAF 
carry with them substantial kinetic energy. This en- 
ergy will be dumped into the external medium and 
will have important consequences. In the context of 
galaxy formation, there has been discussion recently 
of the so-called "radio mode" of accretion (Croton 
ct al. 2006) in which outflowing energy from an ac- 
creting SMBH in the galactic nucleus heats up the 
surrounding medium. This kind of AGN feedback 
can lead to various effects such as shutting off accre- 
tion and/or star formation (Di Matteo et al. 2005; 
Hopkins et al. 2006a). 

It is worth noting that this radio mode is nothing 
other than the ADAF mode of accretion reviewed in 
this article. Investigators of AGN feedback may find 
it profitable to study the considerable work that has 
been done on ADAFs over the last fifteen years. 

4. The Black Hole Event Horizon 

The first BH, Cygnus X-l, was identified and es- 
tablished in 1972 via a measurement of its mass, 
which was shown to be too large for a neutron star 
(NS) . The surest evidence for the existence of BHs 
continues to be through dynamical mass measure- 
ments. We now know of 20 additional compact bi- 
nary X-ray sources (McClintock & Rcmillard 2006; 
Orosz et al. 2007) with primaries that are too mas- 
sive to be a NS or any stable assembly of cold degen- 
erate matter, assuming that GR is valid. Similarly, 
dynamical data have established the existence of su- 
permassive BHs, most notably in the nucleus of our 
Milky Way Galacy (Schodel et al. 2002: Ghcz et al. 
2005a) and in NGC 4258 (Miyoshi et al. 1995). 

Are these compact objects genuine BHs - pock- 
ets of fully collapsed matter that are walled off from 
sight by self gravity and that, like a shadow, reveal 
no detail — or are they exotic objects that have no 
event horizons but manage to masquerade as BHs? 
Most astrophysicists believe that they are genuine 
BHs. There are several reasons for this confidence, 
the most important being that BHs are an almost 
inevitable prediction of GR. However, this argument 
is circular because it presumes that GR is the correct 
theory of gravity. Furthermore, it ignores the many 
exotic alternatives to BHs that have been suggested. 
Thus, the collapsed objects that we refer to through- 
out as "black holes" are strictly speaking dynam- 
ical BH candidates. The current evidence for BHs 
is not decisive, nor can dynamical measurements be 
expected to make it so. We now consider some ap- 



proaches aimed at establishing that these dynamical 
BHs are genuine. 

The defining property of a BH is its event horizon. 
Demonstrating the existence of this immaterial sur- 
face would be the certain way to prove the reality of 
BHs. Unfortunately, unlike any ordinary astronom- 
ical body such as a planet or a star of any kind, it is 
quite impossible to detect radiation from the event 
horizon's surface of infinite redshift. (Hawking ra- 
diation is negligibly weak for massive astrophysical 
BHs.) Nevertheless, despite the complete absence 
of any emitted radiation, it is possible to marshal 
strong circumstantial evidence for the reality of the 
event horizon. The fruitful approaches described be- 
low are based on comparing X-ray binary systems 
that contain BH primaries with very similar systems 
that contain NS primaries. Such investigations are 
motivated by the simple fact that the termination 
of an accretion flow at the hard surface of a NS has 
observational consequences. 

The earliest and strongest evidence for the event 
horizon is based on the faintness in quiescence of 
BH transient systems relative to comparable NS sys- 
tems. In §4. 1 , we discuss the physical arguments that 
underpin this evidence, and in §4.2 we present the 
comparative luminosity data for BHs and NSs. Al- 
ternative explanations for the lower luminosities of 
the BH systems, which do not involve the event hori- 
zon, are considered in §4.3. In §4.4, we present three 
independent and additional arguments for the exis- 
tence of the event horizon, which are again rooted 
in comparing BH and NS X-ray binaries. In §4.5, we 
discuss the extreme faintness and properties of Sgr 
A* and the evidence that this supermassive BH has 
an event horizon, and in §4.6 we argue that even very 
exotic objects with enormously strong surface grav- 
ity (e.g., gravastars) cannot escape our arguments 
for the event horizon. 

4.1. Accretion and the Event Horizon 

A test particle in a circular orbit at radius R 
around a mass M has, in the Newtonian limit, ki- 
netic energy per unit mass equal to GM/2R, and 
potential energy equal to —GM/ R. In the context of 
a gaseous accretion disk, this means that a gas blob 
at radius R retains 50% of its potential energy as 
kinetic energy The remaining 50% was transformed 
into thermal energy during the viscous accretion of 
the blob to its current radius (Frank et al. 2002). 

We now turn to consider accretion on to the ma- 
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terial surface of a NS. If the mass M at the center 
of an accretion disk has a surface, then the kinetic 
energy in the accreting gas will be converted to ther- 
mal energy in a viscous boundary layer and radiated 
(Frank et al. 2002). In addition, the residual ther- 
mal energy that the gas possesses when it reaches 
the inner edge of the disk will also be radiated from 
the surface. 

If accretion occurs via an ADAF, the luminosity 
from the accretion disk L acc and that from the stellar 
surface L SUI f will satisfy 

L acc <GMM/2R*, (14) 
L surf « GMM/R* » L acc , (15) 

where i?* is the radius of the stellar surface. The ac- 
cretion luminosity is small because the flow is radia- 
tively inefficient. Therefore, essentially all the po- 
tential energy of the accreting gas remains in the 
gas in the form of thermal and kinetic energy. If the 
central object has a surface, e.g., it is a NS, the total 
luminosity we observe will be equal to L sm f. How- 
ever, if the object has an event horizon, i.e., it is a 
BH, there will be no radiation from a stellar surface 
and the luminosity we observe will only be equal to 
L acc . Thus, we expect 

ADAF : L NS » L BH . (16) 

Figure 7 shows schematically the luminosity dif- 
ference we predict between a NS and a BH. This plot 
is based on the accretion efficiency estimate shown 
in Fig. 4. Especially at low mass accretion rates, say 
M < 10~ 4 MEdd (deep quiescent state), we expect 
a huge luminosity difference between NSs and BHs. 
Efforts to test this prediction are described in §4.2. 

If the accretion disk is radiatively efficient, then 
the gas has little thermal energy when it reaches the 
inner edge of the disk. In this case, we expectOO 

L acc nGMM/2R*, (17) 
L surf « GMM/2R* ~ L acc . (18) 

The above result is based on a Newtonian analysis 
and is okay so long as the central object has a radius 



9 The result given here is for a non-spinning accretor. If 
the central mass is spinning, the boundary layer luminos- 
ity from the surface will be somewhat smaller. However, it 
requires very fine tuning to make L auI { <g L acc . Moreover, 
if the central mass spins too close to the "break-up" limit, 
the extra luminosity due to the surface can actually exceed 
GMM/2R, by a large factor (Popham & Narayan 1991). 
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Fig. 7. Predicted luminosities of accreting BHs and NSs as 
a function of mass accretion rate. This plot is based on 
the accretion efficiency prescription shown in Fig. 4. Be- 
cause a NS has a surface, the total luminosity is always 
~ GM N gM/i? N g ~ O.lMe 2 , regardless of whether the ac- 
cretion flow is radiatively efficient or not. In the case of a 
BH, however, the radiative efficiency r\ of the accretion flow 
is of paramount importance, and the observed luminosity 
becomes -C O.lMe 2 when the accretion rate is highly sub- 
-Eddington. 

R* ~ -Risco, the radius of the ISCO. For more com- 
pact objects, we must allow for the fact that, inside 
the ISCO, the accreting gas no longer spirals in by 
viscosity but free-falls in the gravitational potential 
of the central mass. Now we have 

L acc «0.1Mc 2 , (19) 
L surf « GMM /R» — L acc > L acc . (20) 

The accretion luminosity is limited by the binding 
energy of the gas at the ISCO, which gives a radia- 
tive efficiency 77 ~ 0.1. On the other hand, the to- 
tal energy budget is GMM/R*, which means that 
a larger fraction of the luminosity is released at 
the surface. (The free-falling gas inside the ISCO, 
crashes on the surface and releases its energy in a 
shock.) In the limit when R* is arbitrarily close to 
Rs (cf., the discussion of gravastars in §4.6), the to- 
tal luminosity is equal to Mc 2 , i.e., the entire rest 
mass energy of the accreting gas is converted to ra- 
diation (100% radiative efficiency). Note that all lu- 
minosities discussed here refer to measurements by 
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an observer at infinity. 

Combining the results in equations (fl"5]) , (fT8|) and 
(|20| . we see that, regardless of the radius of the ac- 
cretor and whether the accretion flow is radiatively 
efficient or inefficient, we expect 

Any Accretion Flow : L SUI { > L acc . (21) 

Although this relation is weaker than (fT"5|) . it can still 
be used in favorable cases to test for the presence of 
an event horizon (§4.5). 

4.2. Evidence for the Event Horizon: Luminosities 
of Quiescent X-ray Binaries 

BH and NS transient binary systems are very sim- 
ilar in many respects, and it is reasonable to ex- 
pect that their mass accretion rates and luminosi- 
ties would be comparable under similar conditions. 
There is, however, one important qualitative differ- 
ence between the two kinds of object — NSs have 
surfaces whereas BHs do not. Often, this difference 
is not important. However, as discussed in §4.1 (see 
eg. fT6|) . when accretion occurs via an ADAF, a NS 
binary should be much more luminous than a BH 
binary (NY95b). The difference will be especially 
large in the quiescent state, when the accretion flow 
is radiatively extremely inefficient (Fig. 7). 

Narayan, Garcia & McClintock (1997b) and Gar- 
cia, McClintock & Narayan (1998) collected avail- 
able X-ray data on quiescent NS and BH transients 
and showed that BH systems are consistently fainter 
than NS systems. This was the first indication that 
the objects that astronomers call "black holes" are 
indeed genuine BHs with event horizons. 

Lasota & Hameury (1998) made the important 
point that it is necessary to compare quiescent NS 
and BH transients at similar mass accretion rates, 
and the surest way to achieve this is to plot lumi- 
nosities as a function of the binary orbital period 
-Porb (the relevant arguments are outlined in Menou 
et al. 1999b; Lasota 2000, 2008). Since 1999, this has 
been the standard way of plotting the data (Menou 
et al. 1999a; Garcia et al. 2001; Narayan, Garcia & 
McClintock 2002a; Hameury et al. 2003; McClintock 
et al. 2004). Figure 8 shows the current status of 
the comparison, with the luminosities of quiescent 
BH and NS transients plotted in (most-appropriate) 
Eddington units. Fig. 9 shows the same data with- 
out the luminosities being scaled. It is clear that, for 
comparable orbital periods, the BH systems are 2 to 
3 orders of magnitude fainter than their NS cousins. 
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Fig. 8. Eddington-scaled luminosities (0.5-10 keV) of BH 
transients (filled circles) and NS transients (open circles) 
versus the orbital period (see Garcia et al. 2001; McClintock 
et al. 2004). Only the lowest quiescent detections or Chan- 
dra/XMM upper limits are shown. The plot shows all sys- 
tems with known orbital periods, which have optical coun- 
terparts and good distance estimates. The diagonal hatched 
areas delineate the regions occupied by the two classes of 
sources and indicate the observed dependence of luminosity 
on orbital period. Note that the BH systems are on average 
nearly 3 orders of magnitude fainter than the NS systems 
with similar orbital periods. 

As discussed in §4.2, a radiatively inefficient 
ADAF provides a natural explanation for the large 
luminosity deficit of the BH systems. In this model, 
the gas approaches the center with a large amount 
of thermal energy. A BH is dim because the bulk of 
this thermal energy is trapped in the advective flow, 
passes through the event horizon, and is lost from 
sight. On the other hand, a NS is bright because 
the thermal energy is radiated from its surface. 

4.3. Other Explanations for the Relative Faintness 
of Quiescent Black Hole Binaries 

A number of alternative models have been put 
forward in an attempt to rationalize the luminosity 
differences between quiescent NS and BH systems. 
Notice below that some of these models generate X- 
ray emission without invoking accretion at all. We 
first consider a recent challenge centered on the dis- 
covery of the extremely low-luminosity NS transient 
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Fig. 9. Same as Fig. 8 except that shown here are the ob- 
served luminosities without Eddington scaling. In this rep- 
resentation, the BH systems are on average a factor of ~ 100 
fainter than NS systems with similar orbital periods. 

1H 1905+00. Next, we discuss the possibility that 
the bulk of the accretion power is channeled into a 
steady jet rather than into X-ray emission. We then 
consider several diverse models that are discussed in 
further detail in Narayan et al. (2002a). 

The case of IE 1905+00: Jonker et al. (2007) 
claim that the extremely low X-ray luminosity of the 
NS transient and type I burst source 1H 1905+00 
(hereafter, H1905) undermines the evidence summa- 
rized in §4.2 for the existence of the event horizon. 
The luminosity of this NS (w 2 x 10 30 ergs" 1 ; D = 
10 kpc), is the same as the luminosity of A0620-00 
and several other short-period (P or b ~ 4— 12 hr) BH 
systems. Based on this result, Jonker et al. assert 
that the evidence for event horizons is "unproven." 
However, they ignore a key point of our argument: 
Namely, as discussed above and illustrated in Fig- 
ures 8 and 9, the case for event horizons depends 
critically on comparing BH and NS systems with 
similar orbital periods. 

We first note that the orbital period of H1905 is 
unknown. More importantly, this unknown period 
is believed to be very short - so short that there 
are no BH systems with comparable periods. This 
conclusion is based on deep optical imaging data 
and the lack of a counterpart. Jonker et al. (2007) 
conclude that the secondary "can only be a brown 



or a white dwarf," and that the system is probably 
an ultracompact binary. Hence the orbital period is 
likely to be tens of minutes, far less than the shortest 
BH binary period of 4.1 hr. 

Thus, there is no BH system comparable to 
H1905. And there is no way to usefully predict the 
mass accretion rate, as Jonker et al. conclude in the 
paper's final sentence. The rough trend of luminos- 
ity with period indicated in Figure 9 for NSs might 
imply a very low luminosity for H1905, as observed, 
but such an extrapolation is unwarranted. See La- 
sota (2007, 2008) for a more detailed discussion of 
this and other issues. 

Jonker et al. (2007) additionally argue that an 
unknown amount of mass transfered from the sec- 
ondary might be lost from the system in winds or 
jets. Their one example of a wind, viz., an extraor- 
dinary mass-loss episode observed in GRO J1655-40 
(Miller et al. 2006c) when the source was in out- 
burst and had a luminosity ~ 10 6 times the quies- 
cent level, is quite inappropriate. Also, their men- 
tion of possible mass expulsion for quiescent NSs via 
the propeller mechanism (Lasota & Hameury 1998; 
Menou et al. 1999a), if effective, would reduce the 
luminosities of NSs and would only strengthen our 
argument. 

Jet outflows: Fender et al. (2003) argue that tran- 
sient BHs at low accretion rates (L/L^dd < 10~ 4 ) 
"should enter 'jet-dominated' states," in which the 
majority of the accretion power drives a radiativcly- 
inefficicnt jet. As we have seen in §3.6, there is good 
evidence for this. For instance, the presence of a ra- 
dio jet has been reasonably well established in qui- 
escence for A0620-00 (Gallo et al. 2006), the closest 
and one of the least luminous of the BHs in question. 
Fender et al. (2003) appeal to the empirical result 
that BHs are ~ 30 times as 'radio loud' as NSs at 
similar accretion rates to deduce that quiescent BHs 
should be ~ 100 times less luminous in X-rays than 
quiescent NSs. Their argument is a little convoluted 
since, if the radio and X-ray emission in quiescence 
come from the jet (see §3.6), one would think the 
BHs would be <~ 30 times more luminous, not 100 
times less luminous. Therefore, naively, the jet ar- 
gument only makes the discrepancy in Figs. 8 and 
9 more severe. Fender et al. (2003) get around this 
difficulty by postulating, in addition to different jet 
efficiencies, also different origins for the X-rays seen 
in quiescent BHs and NSs. 

It is unsatisfying that the reason for the difference 
in jet activity between BHs and NSs "remains un- 
clear." At bottom, the authors are suggesting that in 
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quiescence the BHs are in the 'jet-dominated' regime 
and the NSs are, "if not jet-dominated, close to the 
transition to this regime." A key uncertainty in this 
suggestion is whether the scaling relation that has 
been established for BHs between radio luminosity 
and jet luminosity (X r adio oc L° ,7 ) is also valid for 
NSs, as is assumed. In addition, Occam's razor sug- 
gests that the relative faintness of BHs is unlikely 
to be attributable to their stronger jets: In §4.5, we 
show that Sagittarius A* , which is radiating in quies- 
cence at the same level of Eddington-scalcd luminos- 
ity as A0620-00, does not possess an energetically- 
important jet. 

Coronal emission from BH secondaries: The qui- 
escent X-ray luminosity of the BH systems has been 
attributed to a rotationally-enhanced stellar corona 
in the secondary star (Bildsten and Rutlcdgc 2000, 
but see Lasota 2000). However, as discussed in detail 
by Narayan et al. (2002a), the luminosities of three 
of the systems plotted in Figures 8 and 9 exceed by 
a factor of 6-60 the maximum predicted luminos- 
ity of the coronal model; likewise the three systems 
with adequate data quality show X-ray spectra that 
are harder/hotter than that typically seen in stel- 
lar coronae. Finally, if stellar coronae do contribute 
at some level, then the accretion luminosities of the 
BHs are even lower than our estimates, which would 
further strengthen the evidence for event horizons. 

Incandescent neutron stars: The quiescent lumi- 
nosity of NS transients has been attributed to heat- 
ing of the star's crust during outburst followed by 
cooling in quiescence. This model likewise has prob- 
lems. (For details and references, see Narayan et al. 
2002a.) Briefly, the rapid variability of the prototyp- 
ical NS transient system Cen X-4 is not expected 
in a cooling model and implies that no more than 
about a third of the quiescent luminosity is due to 
crustal cooling. Furthermore, power-law tails that 
carry about half the total luminosity are observed 
for many NS transients (e.g., Cen X-4 and Aql X- 
1). These are unlikely to arise from a cooling NS 
surface, but could easily be produced by accretion. 
Finally, strong evidence for continued accretion in 
quiescence comes from optical variability, which is 
widely observed for these quiescent systems. 

Pulsar wind/shock emission: In another accretion- 
less model, the NS transient switches to a radio 
pulsar-like mode in quiescence (Campana & Stella 
2000). The X-ray luminosity is expected to be of or- 
der the "pulsar shock" luminosity <~ 2 x 10 32 erg s" 1 , 
which is close to the observed level. Both the ob- 
served power-law and thermal components of emis- 



sion (see the previous paragraph) are naturally ex- 
plained by this model: the former component is pro- 
duced by the shock and the latter is radiated from 
the NS surface. The lack of any significant periodic- 
ity in the quiescent emission (in any electromagnetic 
band) could be a problem for this model. Millisecond 
X-ray pulsars such as SAX J1808. 4-3658 do show 
periodicities in outburst, but that emission is clearly 
the result of accretion, not a pulsar wind/shock. 

Optical/ 'UV 'luminosity: Campana & Stella (2000) 
note that the comparison shown in Figs. 8 and 9 as- 
sumes the X-ray luminosity is an accurate measure 
of the accretion rate near the BH or NS. They ar- 
gue that the optical and UV luminosity also origi- 
nates near the central object and should therefore 
be included in the comparison. The non-stellar op- 
tical/UV luminosity is much greater than the X-ray 
luminosity. When it is included, the difference be- 
tween the BH and NS systems largely disappears. 
However, as detailed in Narayan et al. (2002a), there 
are some problems with this argument. The level of 
optical/UV emission generated in the inner region 
of ADAFs is strongly suppressed by winds and con- 
vection, as evidenced by a comparison of the lumi- 
nosities of NSs and white dwarfs (Loeb, Narayan & 
Raymond 2001). It thus appears unlikely that the 
optical and UV emission is generated in the hot gas 
close to the accretor (see Shahbaz et al. 2005). In 
the case of CVs, it has been established that a large 
fraction of the optical/UV emission comes from the 
"hot spot," and it is quite reasonable to expect that 
this is true for BH and NS systems as well (there is 
some evidence to support this; Narayan et al. 2002a). 
Further observations are needed to determine the 
origin of optical/UV emission in these systems. 

4.4. Further Evidence for the Event Horizon in 
BHBs 

In addition to the arguments discussed above, we 
briefly summarize three additional lines of evidence 
for the existence of event horizons. All three are 
based on comparisons between BH and NS X-ray bi- 
naries. As in the examples above, the first argument 
considers quiescent systems, whereas the latter two 
consider active states of accretion. 

In quiescence, a soft component of thermal emis- 
sion is very commonly observed from the surfaces 
of accreting NSs, which is widely attributed to ei- 
ther deep crustal heating (Brown, Bildsten & Rut- 
ledge 1998) or to accretion. No such component is 
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present in the spectrum of the quiescent BHB XTE 
J1118+480 (hereafter J1118), as one would expect 
if the compact X-ray source is a bona fide BH that 
possesses an event horizon (McClintock et al. 2004). 
Because of the remarkably low column density to 
J1118 (iV H w 1.2 x 10 20 cm" 2 ) the limit on a hypo- 
thetical thermal source is very strong (feT^ < 0.011 
keV); it is in fact a factor of ~ 25 lower in flux than 
the emission predicted by the theory of deep crustal 
heating, assuming that J1118 has a material surface 
analogous to that of NSs. Likewise, there is no evi- 
dence that accretion is occurring in quiescence onto 
the surface of J1118, which is the mechanism often 
invoked to explain the far greater thermal luminosi- 
ties of NSs. The simplest explanation for the absence 
of any thermal emission is that J1118 lacks a mate- 
rial surface and possesses an event horizon. 

Type I X-ray bursts are very common in NS X-ray 
binaries, but no type I burst has been seen among 
the BH systems. A model developed by Narayan & 
Heyl (2002, 2003), which reproduces the gross obser- 
vational trends of bursts in NS systems, shows that, 
if the dynamical BHs have surfaces, they should ex- 
hibit instabilities similar to those that lead to type 
I bursts on NSs. Remillard et al. (2006), following 
earlier work by Tournear et al. (2003), measured the 
rates of type I X-ray bursts from a sample of 37 non- 
pulsing X-ray transients observed with RXTE dur- 
ing 1996-2004. Among the NS sources, they found 
135 type I bursts in 3.7 Ms of PCA exposures (13 
sources) with a burst rate function consistent with 
the Narayan & Heyl model. However, for the BH 
group (18 sources), they found no confirmed type 
I bursts in 6.5 Ms of exposure. Their upper limit 
on the incidence of burst activity in these sources 
is inconsistent with the model predictions at a high 
statistical significance if the accretors in BHBs have 
solid surfaces. The results provide strong indirect 
evidence for BH event horizons, and it would appear 
that the evidence can be refuted only by invoking 
rather exotic physics. 

Likewise drawing on the extensive archive of 
RXTE data, Done and Gierlinski (2003) have ex- 
amined the patterns of X-ray spectral evolution 
of active BH and NS sources and identified a dis- 
tinct type of soft spectrum that is occasionally 
observed only in the BH sources. They attribute 
this spectrum to thermal emission from the inner 
accretion disk (corresponding to the high state dis- 
cussed earlier). They then argue that NSs with a 
similar accretion rate cannot exhibit such a sim- 
ple, low-temperature spectrum because they would 



have a second component in the emission from the 
boundary layer where accreting matter impacts the 
stellar surface. They present a thermal/nonthermal 
Comptonization model for the boundary layer emis- 
sion which has significant uncertainty, given that 
so many details of the accretion physics are com- 
plex and poorly understood. Nevertheless, Done & 
Gierlinski (2003) do appear to have identified a sys- 
tematic difference in the X-ray spectra of accreting 
BHs and NSs, and it is surely worth pursuing this 
signature in the effort to amass evidence for the 
reality of event horizons. 

4.5. Event Horizon in Sagittarius A* 

The supermassive BH in Sgr A* has a mass 
M ~ 4 X 1O 6 M (Schodel et al. 2002; Ghez et 
al. 2005a) but an accretion luminosity of only 
~ 10 36 ergs -1 ; thus the source is highly sub- 
Eddington, L < 10 _8 LEdd- From our earlier discus- 
sion, the accretion flow must be in the form of an 
ADAF, i.e., the radiating gas must be extremely hot 
and optically thin. This expectation is confirmed 
by the 10 10 K brightness temperature of the ra- 
dio/millimeter emission (Shen et al. 2005) and the 
fact that most of the emission is in this band rather 
than at frequencies v ~ kT/h (X-ray/7-ray band). 

The ADAF model has been successfully applied to 
Sgr A* (Narayan et al. 1995, 1998a; Manmoto et al. 
1997; Mahadevan 1998; Yuan et al. 2003; to mention 
a few). If we assume that the accretion flow is radia- 
tively very inefficient, then it is easy to make a strong 
case for Sgr A* not having a surface (Narayan et 
al. 1998a). However, as Broderick & Narayan (2006, 
2007) showed, we can argue for the presence of an 
event horizon even without assuming an ADAF. 

If Sgr A* does not have an event horizon, but 
has a surface, then any emission from the surface 
will be blackbody-like. This is because we expect 
the surface to be optically thick. However, the 
radio/millimeter emission mentioned above can- 
not be from this surface because of its incredibly 
high 10 10 K brightness temperature. Thus, the ra- 
dio/millimeter radiation is from the accretion flow, 
and its luminosity gives a strict lower bound on 
L aC c- By equation (|21|) then, we expect a luminosity 
of at least L S urf ~ 10 36 ergs -1 from the surface. The 
surface radiation would be thermal and blackbody- 
like and should come out in the infrared (as easily 
shown, given the luminosity and the likely area of 
the surface). However, there is no sign of it! 
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Fig. 10. The four solid lines show independent upper limits 
on the mass accretion rate at the surface of Sgr A* (assuming 
the source has a surface) as a function of the surface radius R. 
Each limiting curve is derived from a limit on the quiescent 
flux of Sgr A* in an infrared band. The hatched area at 
the top labeled "Typical RIAF Range" corresponds to the 
mass accretion rate in typical ADAF models of Sgr A* (e.g., 
Yuan et al. 2003). The horizontal dashed line represents 
the minimum accretion rate needed to power the bolomctric 
luminosity of Sgr A*. (From Broderick & Narayan 2006) 

Figure 10 shows constraints on the mass accretion 
rate in Sgr A* . The horizontal dashed line is the min- 
imum accretion rate M acc in Sgr A* ; this is the rate 
needed, even with a radiatively efficient flow (rj ~ 
0.1), just to power the observed radio/millimeter ra- 
diation. The four solid lines show the maximum mass 
accretion rate M S urf allowed if Sgr A* has a spher- 
ical surface of radius R. Each curve corresponds to 
a measured limit on the steady quiescent infrared 
flux in a particular band (Stolovy et al. 2003; Clenet 
ct al. 2004; Ghez et al. 2005b) and provides an in- 
dependent upper limit on M sur f . For radii less than 
l5GM/c 2 (the upper limit on the size of Sgr A*; 
Shen et al. 2005), we see that all four bands give up- 
per limits on M sur f that are far below the minimum 
M acc ; in fact, the discrepancy is larger than a factor 
of 100 in the 3.8 fim band. Note that the discrepancy 
would be much larger if we assumed that the accre- 
tion flow is radiatively inefficient (corresponding to 
the hatched region in Fig. 10). 

The only way to avoid the large discrepancy illus- 
trated in Fig. 10 is to give up the assumption that 



Sgr A* has a surface. Obviously, if the source has 
an event horizon, then we do not expect any surface 
radiation, and there is no problem. 

What if Sgr A* ejects all the accreting mass via a 
jet before the gas reaches the surface? Is this a viable 
explanation for the lack of surface emission? We can 
easily rule out this possibility. Recall that the source 
of energy in an accretion system is gravity. At a bare 
minimum, we know that matter is being accreted at 
a rate M acc (shown by the horizontal dashed line in 
Fig. 10) in order to produce the observed radiation. 
All of this mass has to fall into the potential well in 
order to release its energy. If Sgr A* has a jet with 
a kinetic luminosity £j et , the energy for this must 
also come from accretion. We will then require an 
even larger M acc , and hence a larger M sur f , since the 
accreting gas now has to power both the radiation 
and the jet. The discrepancy with the observed lim- 
its on the quiescent infrared flux would then be even 
larger. 

4.6. Effects of Very Strong Gravity 

The arguments for the event horizon presented in 
§§4.2, 4.5 are based on Newtonian ideas. Some au- 
thors (e.g., Abramowicz, Kluzniak & Lasota 2002) 
have questioned whether the arguments might be 
substantially modified by strong gravity in the vicin- 
ity of the compact accretor. 

Buchdahl (1959) showed, for a wide class of rea- 
sonable equations of state, that the smallest radius 
allowed for a compact non-rotating object is i? m in = 
(9/8)i?g. An object with this limiting radius has a 
gravitational redshift from its surface z gT — 2. At 
such modest redshifts, we do not expect the effects 
of strong gravity to be particularly large, and so the 
arguments presented in §§4.2, 4.5 will continue to 
hold. 

Recently, however, a new class of solutions has 
been discussed, which goes variously under the name 
of "gravastar" and "dark energy star" (Mazur & 
Mottola 2001; Chapline et al. 2003; Visser & Wilt- 
shire 2004; Carter 2005; Lobo 2006). In this model, 
the radius R of an object of mass M is allowed to be 
arbitrarily close to R s : AR = R - R s < Rs- The 
surface redshift can then be arbitrarily large: 

1 + z gI « (Rs/AR) 1 / 2 > 1. (22) 

Although the gravastar model is very artificia l 10 1 , 
it is nevertheless interesting to ask whether such a 

Broderick & Narayan first submitted their 2007 paper to 
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model, which has no event horizon, can explain the 
observations described in §§4.2, 4.5 (Abramowicz et 
al. 2002). Broderick & Narayan (2006, 2007) show 
that it cannot. We briefly summarize the arguments 
here. 

A large value of z gr means that any radiation 
emitted from the surface is redshifted greatly before 
it reaches the observer at infinity. This looks like 
an easy way of hiding the surface luminosity. How- 
ever, a simple energy conservation argument shows 
otherwise. Assuming steady state, in the limit as 
AR/Rs — > 0, the total luminosity observed at infin- 
ity, L acc + L SU rf , must be equal to Mc 2 (the rate of 
accretion of rest mass energy). Moreover, the effec- 
tive radius of the source as viewed by a distant ob- 
server is (3\/3/2)-Rs = 2.6i?s, and so the surface ra- 
diation will be in the X-ray band for accreting BHBs 
(§4.2) and in the infrared for Sgr A* (§4.5), exactly 
as in the Newtonian analysis. 

Could the large gravitational redshift cause a large 
delay in the signals from the surface of the star, 
and could this be why we do not see the radiation? 
The extra delay due to relativity is easily estimated 
by considering null geodesies in the Schwarzschild 
metric: 

At = (R s /c) \n(R s /AR). (23) 

Since the dependence is only logarithmic, the ex- 
tra delay is not significant. For instance, even if we 
take AR to be comparable to the Planck length ~ 
10 -33 cm (the smallest length we can legitimately 
consider), the delay is only <~ 10 ms for a BHB and 
- 1000 s for Sgr A*. 

We assumed that the radiation from the surface 
would have a blackbody-like spectrum. Is this likely? 
Actually, it is virtually guaranteed as AR/Rs — > 0. 
When the radius of an object is very close to the 
Schwarzschild radius, most rays emitted from the 
surface are bent back on to the surface, and only a 
tiny fraction of rays escapes to infinity. The solid an- 
gle corresponding to escape is Q esc ~ AR/Rs <C 1. 
In this limit, the surface behaves just like a furnace 
with a pinhole (the textbook example of a black- 
body!), and so the escaping radiation is certain to 
have a nearly perfect blackbody spectrum. 



Phys. Rev. Lett., but the paper was not even sent out for 
review because, in the words of the editor, "the gravastar 
model of Mazur and Mottola is not considered by the com- 
munity to be a viable alternative." B&N then submitted 
the paper to Class. Quantum Grav., where most papers on 
gravastars and dark energy stars are published. 



Could the surface luminosity escape in the form 
of particles rather than radiation? By the blackbody 
argument given above, whatever escapes must be in 
thermodynamic equilibrium at the temperature 
observed at infinity. For the sources of interest to us, 
fcToo - eV (Sgr A*) and keV (BHBs). In thermody- 
namic equilibrium, the only particles with any sig- 
nificant number density at these temperatures are 
neutrinos. Even if we include all three species of neu- 
trinos, the radiation flux is reduced by a factor of 
only 8/29 (Broderick & Narayan 2007). This is a 
small correction compared to the large discrepancies 
> 100 that we described in §§4.2, 4.5. 

Finally, we note that our arguments for the event 
horizon are based on a steady state assumption. 
Specifically, we assume that the surface luminosity 
is proportional to the average mass accretion rate 
on the surface. In the gravastar model, in particular, 
one can imagine scenarios in which steady state is 
not reached. However, Broderick & Narayan (2007) 
show that this explanation can be ruled out, at least 
with current gravastar models. 

We thus conclude that strong gravity effects are 
unable to weaken our arguments for the presence of 
an event horizon in quiescent BHBs and in Sgr A*. 
We must look to more conventional astrophysical 
explanations. The only idea that we find somewhat 
plausible is that the radiation we observe from qui- 
escent X-ray binaries and Sgr A* is not produced 
by accretion at all, but by some other process. This 
would undercut all our arguments since we assume, 
as an article of faith, that the radiation we observe 
is powered by gravity through accretion. Given the 
success of the accretion paradigm in explaining a 
vast body of observations on a variety of BH systems 
in many different spectral states, it seems rather ex- 
treme to abandon the idea of accretion in just those 
particular sources where we find evidence for the 
presence of an event horizon. 

5. Conclusion 

The aim of this article is two-fold: To give an 
updated and current account of the ADAF model, 
with an emphasis on applications to accreting BHs 
(§3), and to review the considerable body of evi- 
dence presently available for the presence of event 
horizons in astrophysical BHs (§4). 

Observations of BH binaries (BHBs) and active 
galactic nuclei (AGN) indicate that, at luminosities 
of a few percent or less of Eddington, accretion oc- 
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curs via a very different mode than the standard thin 
accretion disk. At these luminosities, sources have 
hard X-ray spectra, quite unlike the soft blackbody- 
like spectra seen in more luminous sources. Prior to 
the establishment of the ADAF model, observations 
in the hard state were modeled in an ad hoc way 
using empirical thermal Comptonization models. In 
the mid-1990s, the ADAF solution was shown to 
have precisely the densities, temperatures, radiative 
(in) efficiencies - and stability - required to provide 
a physical description of the observations. More- 
over, since the ADAF model is essentially mass- 
independent, observations of BHBs and AGN are 
both explained with more-or-less the same ADAF 
model. The model gives satisfactory results for a 
wide range of luminosity, from about 10 _1 down to 
10~ 8 of Eddington (below which there are no obser- 
vations). 

The defining characteristic of an ADAF is that 
the gas is not radiatively efficient. A significant frac- 
tion of the energy released by viscous dissipation is 
retained in the accreting gas and advected with the 
flow. The trapped thermal energy causes the accret- 
ing gas to be weakly bound to the BH. Based on 
this property, it was predicted already in the earliest 
ADAF papers that sources in the ADAF state would 
have strong winds and jets. Nonthermal radio emis- 
sion has been detected in recent years from BHBs in 
the hard state and quiescent state, as well as from 
all AGN with luminosities below ~ 1% of Edding- 
ton. These sources have spectra consistent with the 
ADAF model, thus confirming a strong connection 
between ADAFs and jets/outflows. A currently ac- 
tive topic of research is the role of AGN feedback 
on galaxy formation. Research in this area can now 
draw on the extensive literature on ADAFs. 

Of special interest in this work are those accret- 
ing objects for which there exist estimates or con- 
straints on the two system specific parameters: BH 
mass M and mass accretion rate M. Specifically, we 
are referring here to Saggitarius A* and a selected 
sample of quiescent BH and NS X-ray binaries. We 
have featured these two examples because they pro- 
vide strong observational evidence for the existence 
of the event horizon. 

In the case of the binaries, one can make a plausi- 
ble argument that, for comparable orbital periods, 
the mass accretion rates and luminosities of both 
types of systems should be comparable. Remarkably, 
however, the BH systems are observed to be dimmer 
by a factor of <~ 100 — 1000. We review a wide range 
of unsatisfactory attempts to explain this large lu- 



minosity difference. In contrast, the ADAF model 
provides an entirely straightforward explanation for 
the faintness of quiescent BHs: A NS must radiate 
the trapped thermal radiation that is advected with 
the accretion flow and rains down on its surface, 
whereas the BH hides the energy behind its event 
horizon. Additional evidence for the event horizon 
is provided in further comparisons of NSs and BHs 
in outburst: The BHs lack type I X-ray bursts, and 
they lack a distinctive boundary-layer component of 
emission. Both of these properties are expected if 
the BHs possess event horizons, but hard to explain 
otherwise. 

The supermassive BH in Sgr A* is extraordinar- 
ily quiescent. The observations strongly support the 
existence of an ADAF, and it is consequently easy 
to provide a compelling argument for an event hori- 
zon. However, for Sgr A*, one can make an even 
stronger argument for the lack of a hypothetical ma- 
terial surface without at all invoking an ADAF. The 
large radio/millimeter accretion luminosity of Sgr 
A*, which has a brightness temperature > 10 10 K, 
obviously cannot be emitted by an optically-thick 
surface, and so it must be radiated from the accre- 
tion flow. This establishes a hard lower limit on the 
accretion luminosity and M . Meanwhile, high an- 
gular resolution radio observations constrain the ra- 
dius of the surface to be < 15i?s- This constraint 
and the lower limit on M predict a near-IR flux of 
thermal, blackbody-like surface emission that is far 
above the observed limits. The obvious explanation 
is that there is no material surface, only an event 
horizon. As an added bonus, it is energetically im- 
possible to explain away the lack of surface emission 
by appealing to mass loss in a jet. 

The case for event horizons in both Sgr A* and 
in the stellar-mass BHs is robust against appeals to 
strong gravity or the leading models of exotic stars. 
First, GR effects will be mild for nearly all conven- 
tional models of degenerate stars, whose radii are re- 
stricted to be > 9/8i?s and whose surface redshifts 
are therefore < 2. Secondly, even for an exotic star 
(a gravastar) with an extraordinary surface rcdshift 
of a million or more, the full accretion luminosity 
from its surface will be delivered as X-rays (in bina- 
ries) or infrared (in Sgr A*) to a distant observer. 
Thus, extreme redshifts have practically no effect on 
the argument for the event horizon. 

To play on Carl Sagan's famous comment, there 
will always be an absence of direct evidence for the 
event horizon, but this surely cannot be taken as 
evidence of its absence in nature. On the contrary, 
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and with a pun in mind, many indicators show that 
the event horizon is an inescapable reality. 

The authors thank Bozena Czerny, Jean-Pierre La- 
sota and Feng Yuan for helpful comments. 
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